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We present a calculation of low energy magnetic states of doubly-closed-shell nu- 
clei. Our results have been obtained within the random phase approximation using 
different nucleon-nucleon interactions, having zero- or finite-range and including a 
possible contribution in the tensor channel. 



I. INTRODUCTION 



In the last thirty years electron-nucleus scattering experiments have produced a large 
amount of high precision data, which impose severe constraints on nuclear models and effec- 
tive theories, such as the Random Phase Approximation (RPA). In particular the description 
of low energy excited states within the RPA is known to be very sensitive to the details of 
the effective nucleon-nucleon (NN) interaction used in the calculations. 

We present here a selection of results from a systematic study of the low energy spectra of 
several doubly-closed-shell nuclei we have made within the RPA theory [1]. We have focused, 
in particular, on the unnatural parity states, which are sensitive to the spin, spin-isospin 
and tensor channels of the residual NN interaction 2|. 

In the first step of our project we have employed a purely phenomenological approach, 
using a single-particle mean-field basis generated by a Woods-Saxon well, and constructing 
phenomenological residual NN interactions which reproduce some selected nuclear states. 
In order to study the sensitivity of our results to the details of the residual interaction, 
we have used four different NN interactions, which have zero- and finite-range and may 
include contributions in the tensor channels. We have found some states which are very 



* Email address: chiara.maieron@le.infn.it 



sensitive to the details of the interaction, and in general we have obtained a satisfactory 
description of the low energy states of the nuclei under consideration. We can thus consider 
our phenomenological approach as an "optimal" RPA approach in terms of comparison with 
the experimental data. 

In our study we have then used the RPA amplitudes X and Y to calculate electron 
scattering response functions. As a further independent extension of this approach, we have 
also considered the computation of neutrino cross sections, for which an example will be 
shown in the following. 

We have then proceeded to the second step of our study, performing RPA calculations 
within a self-consistent approach, where the single-particle mean-field basis is obtained by 
means of a Hartree-Fock calculation which uses the same effective NN interaction used in the 
RPA calculation. In particular we have used the Gogny Dl finite-range interaction 
finding, in this case, remarkable disagreement with the experimental spectra. 



II. FORMALISM 



The inputs required by RPA calculations are a set of single-particle energies and wave 
functions and a residual NN interaction. 

In the purely phenomenological calculations, for the former we have used the single- 
particle basis generated by a Woods-Saxon well, whose parameters are taken from the liter- 
El 

ature [6| . For the latter, we have considered a generic residual interaction written as 

Kff(l,2) = «i(ri3) + v^n^p ^,^) 

+ [v 2 (r 12 ) + v%(r 12 ) p a (r ll r 2 )] n ■ r 2 
+ v 3 (r 12 ) o"i ■ o"2 + v 4 (r 12 ) cti • cr 2 ri • r 2 

+ v 5 (r 12 ) Si 2 (fi 2 ) + v 6 (r 12 ) Si 2 (fi 2 ) Ti • r 2 , (1) 

where cr and r are the usual spin and isospin operators and 5i 2 is the tensor operator. As 
suggested by past phenomenological RPA studies 0, Q] we have included density dependent 
terms in the the central and isospin channels, with 

p(ri,r 2 ) = [p(ri)p(r 2 )] 1/2 , (2) 

and using a = 1 in Eq. ([T]). 




FIG. 1: Effective NN interactions as functions of the relative momentum. The dashed lines repre- 
sent the LMtt interaction, the dotted the FRtt and the solid the Gogny Dl. The central channels 
Vi (i=l-4) of the LMtt and FRtt interactions are identical to those of LM and FR, respectively. 

We have considered four different forms of the NN interaction, parametrizing them ac- 
cording to the following criteria: (i) we have chosen a unique set of parameters for all the 
nuclei under investigation, with the exception of the density dependent terms which are dif- 
ferent for each nucleus; (it) the density dependent terms have been set to reproduce the first 
2 + state in 12 C and the first 3 _ states in the other nuclei; (in) the remaining contributions 
in the central and isospin channels have been chosen to get a reasonable description of the 
centroid energy of the isovector giant dipole resonance; (iv) the spin, spin-isospin and tensor 
channels have been adjusted to describe the low energy (below 8 MeV) magnetic spectrum 
of 208 Pb , with particular attention to the 12~ and 1 + states and, in addition, taking care 
that the energy of the first 4" state of 16 is reproduced reasonably. 

Following a Landau-Migdal approach, we have first considered zero-range interactions 
without and with a tensor-isospin channel contribution (LM and LMtt in the following). 

We have then constructed finite-range interactions (FR and FRtt), by keeping the long- 
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FIG. 2: (Color online) Low energy spectrum of 208 Pb calcu 
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ated with the four interactions used in 



this work and compared with the experimental energies 

range behavior of the Argonne vis potential [8] and substituting its short-range part with 
a sum of Gaussians. We have also used Gaussians to parametrize the density dependent 
terms of the interaction and, for the FRtt case, we have obtained the tensor channel terms 
by multiplying the corresponding terms of the interaction by a correlation function 



obtained in variational calculations 



6] 



The behavior of the various interactions we have considered is shown in figure [T] as a 
function of the relative momentum of the interacting pair. 

In the self-consistent ap pro ach we have obtained the single particle basis by solving 



Hartree-Fock equations 19, 10 1, using the same effective NN interaction used for the RPA 

n n ft 

calculations. We have used the Gogny Dl interaction [3|, 14, [5j], shown in figure Q] by the 
solid lines, which has finite-range components in the central, isospin, spin and spin-isospin 
channels and a zero-range density dependent contribution. We have not included the spin- 
orbit term in the RPA calculations. 

For the self-consistent calculations we have included the contributions of both direct and 
exchange matrix elements of the interaction, whereas for the phenomenological approach we 




FIG. 3: (Color online) Right panels: electron scattering transverse responses for the first 12 
states of 208 Pb , versus the effective momentum transfer. Different residual interactions are used 
as indicated. Left panels: the same for the first 10~ states. Experimental data from [jjj]. 

have considered direct terms only, assuming that the effect of the exchange terms is effectively 
included in the choice of the parameters characterizing the various NN interactions. 

Finally we observe that our calculations have been obtained by discretizing the contin- 
uum. We have checked that our results are stable with respect to the parameters charac- 
terizing the continuum discretization and also with respect to the size of the single-particle 
configuration space used for each nucleus. More details on the role of the continuum dis- 
cretization in self-consistent RPA calculations can be found in [2fl. 



III. RESULTS 

We have performed systematic calculations of the low energy spectra of 12 C , 16 , 40 Ca , 
Ca , 90 Zr and 208 Pb . For each state we have also computed electromagnetic response 



functions, and we have compared them with the available experimental data. 

As an example of low energy spectrum, in figure [2] we show the case of 208 Pb . Here we 
present only those experimental states we have been able to identify with those obtained in 
our calculations. We observe that the general agreement is quite good and that, except for 
some cases, the various energies have little sensitivity to the details of the interactions. 

The transverse response functions of the 12~ states of 208 Pb are shown in the right 
panels of figure El We can see that for the higher state at 7.08 MeV (lower right panel) the 
experimental data are rather well reproduced with all NN interactions, which do not produce 
significant differences in the curves. On the other hand, for the lower energy state (upper 
right panel) we observe a very strong dependence on the residual interaction, both when 
finite-range and when tensor channel contributions are included. We remark that a better 
description of this state alone could be obtained with a different choice of the parameters of 
the residual interaction, but this would worsen the global descritpion of the various magnetic 
spectra we have considered. 

In the left panels of the same figure the transverse responses of the 10~ states are also 
shown. It is interesting to notice that, in this case, only the curves which include tensor 
contributions are able to reproduce the second peak shown by the data, for both states. 

An example of the results we have obtained within the self-consistent approach is given in 
table U where the low energy magnetic spectrum of 12 C is presented, and in figure HI where 
we show the corresponding transverse responses of the 1 + isospin-doublet of states. For the 
sake of comparison between the self-consistent results and the purely phenomenological ones, 
we have taken as reference the FRtt case, considered to be the most complete interaction 

TABLE I: Energies of the low- lying magnetic states of 12 C (in MeV). Experimental energies from 
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FIG. 4: Electron scattering transverse responses of the first 1 + states of 12 C , versus the ef- 
fective momentum transfer. The self-consistent results obtained with the Gogny Dl interaction 
are compared with those corresponding to the FRtt phenomenological case. Experimental data 
from QQ. 

we have obtained. 

We see that the magnetic states of 12 C are reproduced rather badly by the self-consistent 
calculations: the response functions agree in shape and magnitude with the data, but the 
order of the states forming isospin doublets is inverted. 

This happens also for 16 , as illustrated, for the case of the 4~ states, in figure El where 
the values of the excitation energies obtained with the Dl and FRtt interactions are also 
reported. 

We have systematically obtained this kind of inversion for all magnetic states of all the 
nuclei we have studied . This result indicates the inadequacy of the Dl interaction in 
isospin-dependent channels, and it may be considered the most important outcome of the 
self-consistent part of our study. 
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FIG. 5: Same as figured] for the first 4 states of 16 . Experimental data from 



IV. EXTENSION OF THE PHENOMENOLOGICAL APPROACH TO 

NEUTRINO SCATTERING 



After obtaining an "optimal" phenomenological RPA approach, by tuning NN interactions 
to reproduce data from electron scattering at best, we have also considered the possibility of 
applying it to the study of low energy neutrino scattering cross-sections, with the purpose 
of studying their sensitivity to the tensor components of the interaction. 

In figure [6] we show the cross sections for the neutral-current neutrino and antineutrino 
excitation of the low energy 0~ states of 16 , as a function of the incident neutrino energy. 
The values of the excitation energies u) are indicated inside the figure. 

We observe that, in this test-case, extremely large differences in the cross sections are 
obtained when tensor channel contributions are included in the effective NN interaction. 
This indicates that the role of the residual interaction in neutrino scattering cross section is 
a very interesting topic, for example in connection with the problem of nuclear uncertainties 
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FIG. 6: Cross sections for the neutral current neutrino and antineutrino excitation of states of 
16 O at excitation energy u, versus the incident (anti)neutrino energy. 



in the detection of supernova neutrinos 



16| 



V. SUMMARY AND CONCLUSIONS 



We have made a systematic study of the low-lying (magnetic) spectra of doubly-closed- 
shell nuclei. Our study indicates that a simultaneous description of all states imposes strong 
constraints on the residual NN interaction. Within a purely phenomenological approach it 
is possible to get a good description of the spectra and of the response functions of most 
of the states, thus obtaining an "optimal" RPA approach. Some states which are not well 
described exhibit a strong sensitivity to some details of the residual interaction, and a deeper 
investigation could be used to obtain further constraints on it. 

On the contrary, the self-consistent calculations we have performed with the Gogny Dl 
interaction produce results in disagreement with the experimental spectra, indicating the 



presence of problems in particular in the isospin dependent parts of the interaction. A 
detailed study can be found in [2j]. 

Finally we have considered an example of a potentially important application of our phe- 
nomenological approach to the calculation of neutrino scattering cross sections, to investigate 
the sensitivity of the latter to the tensor components of the residual interaction. 
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